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Abstract. From the quantum mechanical point of view, the electronic characteristics of quasicrystals are 
determined by the nature of their eigenstates. A practicable way to obtain information about the properties 
of these wave functions is studying the scaling behavior of the generalized inverse participation numbers 
Z q ~ j^~ D i(i~ 1 ) with the system size N. In particular, we investigate d-dimensional quasiperiodic models 
based on different metallic-mean quasiperiodic sequences. We obtain the eigenstates of the one-dimensional 
metallic-mean chains by numerical calculations for a tight-binding model. Higher dimensional solutions of 
the associated generalized labyrinth tiling are then constructed by a product approach from the one- 
dimensional solutions. Numerical results suggest that the relation Dq d = dD^ d holds for these models. 
Using the product structure of the labyrinth tiling we prove that this relation is always satisfied for the 
silver-mean model and that the scaling exponents approach this relation for large system sizes also for the 
other metallic-mean systems. 
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1 Introduction 

Since the discovery of quasicrystals [UEIE] various experi- 
mental investigations have revealed the rather uncommon 
electrical, magnetic and optical characteristics of real qua- 
sicrystals [HE] and motivated extensive research to ob- 
tain a better theoretical understanding of their structure 
and properties. Today many exact results are known for 
one-dimensional quasicrystals [BJ, but the characteristics 
of wave functions in two or three dimensions are under- 
stood to much lesser degree. Also numerical investigations 
often deal only with finite, relatively small systems or pe- 
riodic approximants, which allows only the estimation of 
a few properties of ideal quasicrystals [TBI] . 

An alternative approach is based on the study of d- 
dimcnsional models with separable Hamiltonians such that 
higher dimensional, quantum mechanical solutions can be 
directly derived from the one-dimensional case, allowing 
the numerical consideration of large systems up to 10 8 
sites. The best studied model within this approach is based 
on the Fibonacci sequence [TUllTTirT^] , whereas in this pa- 
per we concentrate on the more general class of metallic- 
mean quasiperiodic sequences [TSirHIITB] . which has been 
less investigated so far. The sequences describe the weak 
and strong couplings of atoms in a quasiperiodic chain and 
the one-dimensional eigenstates are obtained by numeri- 
cal diagonalization of the respective tight-binding Hamil- 
tonian. Higher-dimensional quasiperiodic tilings, the so- 
called generalized labyrinth tilings, are constructed as a 
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direct product of these chains and their eigenenergies and 
eigenstates can be directly calculated by multiplying the 
energies E and wave functions \P of the quasiperiodic chain, 
respectively [15] . 

While for an electron in a one-dimensional quasiperi- 
odic system many examples lead to purely singular con- 
tinuous energy spectra and multifractal eigenstates 17, 18, 
IT5] . in the two-dimensional labyrinth tiling the Lebesgue 
measure of the energy spectrum shows a transition from 
zero to a finite value with increasing coupling parameter v 
[20] . Regardless of this transition of the energy spectrum, 
the numerical results show that the wave functions remain 
multifractals in two and three dimensions and do not be- 
come extended with increasing values of v. To address 
this, we investigate the scaling behavior of the general- 
ized inverse participation numbers Z q of the eigenstates 
and their dependency on the dimension d in more detail. 
Further, we mathematically prove that the generalized di- 
mensions associated to the inverse participation numbers 
Z q satisfy the relation D^ d = AD* d for the silver-mean 
model and asymptotically approach this relation for the 
other metallic-mean systems. 

This paper is organized as follows: In Sec. [2] we intro- 
duce the inflation rule for the quasiperiodic metallic-mean 
chains and discuss the properties of the corresponding one- 
dimensional eigenstates as well as the construction of the 
eigenstates of the associated higher-dimensional labyrinth 
tilings. Section [3] then focuses on the inverse participation 
numbers by discussing numerical results and proving the 
above mentioned relation between the generalized scaling 



2 



Stefanie Thiem and Michael Schreiber: Generalized Inverse Participation Numbers in Metallic-Mean Systems 



exponents D q in different dimensions. This is followed by 
a brief summary of our results. 



2 Eigenstates and Wave Functions of 
Generalized Labyrinth Tilings 

The construction of the labyrinth tiling is based on the so- 
called metallic- mean quasiperiodic sequences, which for a 
parameter 6 are defined by the inflation rule 



V 



„b-l 



(1) 



starting with the symbol w. After a iterations we obtain 
the ath order approximant C a of the quasiperiodic chain. 
The length f a of an approximant C a is also given by the 
recursive rule f a = bf a -i + f a -2 with / = fx = 1. 
Thereby, the ratio of the lengths of two successive iter- 
ants in the limit linia^oo f a /f a -i = A equals different 
metallic means depending on the parameter b with a con- 
tinued fraction representation A = [6] = [6,6,6,...]. This 
leads to the well known Fibonacci sequence for 6 = 1 
with the golden mean Aa u = [1] = (1 + v5)/2, while 
6 = 2 results in the octonacci sequence with silver mean 
AA g = [2] = 1 + \/2 and 6 = 3 corresponds to the bronze 
mean A Bz = [3] = (3 + Vl3)/2 [21]. 



2.1 Eigenstates of metallic-mean chains 

Solving the time-independent Schrodinger equation for the 
quasiperiodic systems with zero on-site potentials 



E l Vi = t l _ 1 M_ 1 +tu+ifi+i , ( 2 ) 



we obtain the discrete energy values E l and the wave func- 
tions I"? -1 ) = represented in the orthogonal 
basis states \l) associated to a vertex The hopping pa- 
rameter t in the Schrodinger equation is given according 
to the quasiperiodic sequence C a with t s = 1 for a strong 
bond and t w = v for a weak bond (0 < v < 1) [711X5] . 
Applying free boundary conditions the number of vertices 
is N a = f a + 1. 

The results show that the eigenvalues are symmetric 
with respect to 0. Thus, for even system sizes N a all energy 
values E have a symmetric value —E, but for odd N a there 
is one state E M — 0, which has no corresponding state. 
Here M — \N a /2\ + 1, assuming that the eigenstates are 
labeled according to increasing eigenenergies. 

Additionally, the eigenfunctions possess a symmetry. 
The eigenstate tf^ with the eigenvalue E and the eigenstate 
^ with the corresponding eigenvalue — E only differ by an 
alternating sign depending on the vertex I according to 



*i = (-1)% 



(3) 



Further, for odd N a = 2M — 1 the eigenvector asso- 



namely W t 

i (cp. na 



M 



vanishes either on all odd or on all even sites 
i.e. 



\Ir M — n 

w l mod 2=0 — U 



and 



i& M+ - n 

W l mod 2=1 — U 



(4) 



2.2 Eigenstates of the generalized labyrinth tilings 



The generalized labyrinth tilings in d dimensions are con- 
structed from d quasiperiodic chains C a perpendicular to 
each other, where the bonds are given by the diagonals. 
Hence, depending on the starting point the grid decom- 
poses into 2 d ~ 1 separate grids. Each of these grids corre- 
sponds to a finite ath order approximant L a of the gen- 
eralized labyrinth tiling C [TBIITS] . In the case of the oc- 
tonacci sequence all 2 d ~ 1 grids are identical. For other in- 
flation rules the different approximants are slightly shifted 
against each other and mainly differ at the boundaries. 
These similarities originate from hidden mirror symme- 
tries in the quasiperiodic chains. While the octonacci chain 
is perfectly mirror symmetric, for the other chains mirror 
symmetry can be achieved by neglecting a few symbols at 
one end of the sequence (e.g. for the Fibonacci chain the 
last two symbols on the right) and by interchanging two 
consecutive bonds for the case 6 > 2 (i.e. a phason flip) 

m 

The eigenstates of the generalized labyrinth tiling in 
d dimensions are then constructed as the product of the 
eigenstates of d one-dimensional chains. For instance in 
two dimensions we use E ij = E ll E 2 ' and ~ 

LTD, L IlL 

[15) . The indices i and j enumerate the eigenvalues E in 
ascending order and I and m represent the coordinates 
of the vertices. However, due to the symmetries of the 
eigenfunctions of Eq. some of the wave functions $\ 3 m 
and the related eigenvalues E l i are identical and then only 
one of them is allowed to be considered. 

In two dimensions for even chain lengths N a the system 
size is N%/2 and a valid combination of eigenstates is 

{^|l<t<JV„/2Al<j<JV } . (5) 
The normalization of the wave functions results in 



lm 



I 771 



(6) 



For odd chain length N a the selection of the eigenstates 
is more complicated due to the special structure of the 
eigenfunction for E M = 0. This results in the set of two- 
dimensional energy values 

{E l] | (1 < i < M A 1 < j < N a ) V 

(i = M A 1 < j < M)} . (7) 

Further, the normalization of the wave functions is given 

by 

fmL i=j=M 



yy/2^l^ m otherwise 



(8) 



The special treatment of the case i = j = M is due to the 
fact that the wave function <Pf^ lives only on one of the 



ciated to the eigenvalue E — has a special structure, two different labyrinth tilings, C or £*. 
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We analogously compute the eigenenergies and wave 
functions in higher dimensions. For even N a the eigen- 
values E l —i k — E % .. .E^E k of the generalized labyrinth 
tiling £f d in d dimensions are described by the set 

{E l - jk \l<i< N a /2 A ... A 1 < j < N a /2A 

l<k<N a }. (9) 

This equation results from the symmetry of the wave func- 
tions given in Eq. ^ . Due to the special structure of the 
eigenfunctions with the eigenvalue E M — 0, for odd N a 
the allowed set of eigenstates is much more complex. For 
instance, in three dimensions a valid combination is given 
by 

{E ljk | (1 < i < M A 1 < j < M A 1 < k < N a ) V 

(i = M A 1 < j < M A 1 < k < M) V (10) 
(1 < i < M A j = M A 1 < k < M)} . 

Further, the normalization constants of the wave functions 
^imn f° r even system sizes N a are 

k = V2 d - 1 V}& : > <P k (11) 

in analogy to two dimensions. For odd N a and three di- 
mensions we obtain 

[vmvi i=j = k = M 

= \ two indices i, j, k equal M (12) 

[2^i<P^ k otherwise. 

For larger approximants the wave functions approach 
each other because the influence of the boundaries and of 
the phason flip for b > 2 vanishes [TB] . This characteristic 
is useful for analytical considerations of properties of the 
eigenstates of the general labyrinth tiling. 

3 Scaling Behavior of Generalized Inverse 
Participation Numbers 

In a d-dimensional system with discrete positions r the 
generalized inverse participation number of an eigenstate 
<P r is defined for a parameter q as 

Z q (4>) = J2\$r\ 29 ■ (13) 

r 

This is associated to the scaling of the gth moment of the 
eigenstates ^ r according to a power law 

Z q (<P) ~ jy-Diti- 1 ) (14) 

with the linear dimension of the system and the gener- 
alized dimensions D q . 

In particular, the inverse of Zi equals the total num- 
ber of sites for which the probability measure of the wave 
function given by |^ r | 2 is significantly different from zero 
and hence is a measure for the localization of the wave 



function <P r . Investigating the scaling behavior of the par- 
ticipation ratio 

p(<p) = 7,2 Jg) ~ N-^ (15) 

in the volume V, the wave function of a localized state is 
characterized by the scaling exponent 7 = 1, and 7 = 
corresponds to an extended state. Intermediate values of 7 
(0 < 7 < 1) indicate fractal eigenstates, which are neither 
extended over the whole system nor completely localized 
at a certain site and show self-similar patterns jT? lf2"2"ll2"5] . 

Within a small energy interval the participation ra- 
tios usually fluctuate over a certain range and the average 
participation ratio is often also energy dependent. For in- 
stance, for octagonal tiling models, for the Penrose tiling 
and for the two- and three-dimensional Rauzy tilings the 
participation ratios are found to be smaller at the cen- 
ter of the energy spectrum than at its edges, while for 
the three-dimensional Ammann-Kramer-Neri tiling the in- 
verse behavior can be observed [7ll8ll2"4 l l2"5 ] . However, to 
get information about the spatial distribution of the wave 
functions for macroscopic systems one has to compute the 
scaling behavior of the participation ratios with the sys- 
tem size. In particular, numerical results for the Penrose 
tiling and the Ammann-Kramer-Neri tiling revealed that 
the scaling exponents "f{E) show a qualitatively similar 
dependence on the energy as the participation ratios [7]. 

In contrast to the former discussed systems, for the 
metallic-mean models the participation ratios and inverse 
participation numbers do not show any significant trend 
over the whole energy spectrum as visualized in Fig. [l] 
This might be related to the same number of nearest 
neighbors for every vertex in these tilings. In order to 
compare the characteristics of the eigenstates for differ- 
ent dimensions we consider the inverse participation num- 
bers {Z2) and the participation ratio (p) averaged over 
the complete energy spectrum. Exemplarily, we show nu- 
merical results for the dependency of the corresponding 
scaling exponents D2 and 7 on the coupling constant v 
for the golden-, silver-, and bronze-mean model in one, 
two and three dimensions in Fig. [2] This indicates that 
both scaling exponents are independent of the dimension 
[16] . The deviations in three dimensions are most likely 
caused by the consideration of relatively small sizes of the 
associated one-dimensional chains with A^, u = 234 and 
A^^ z = 143 sites. Similar results can also be found for the 
scaling exponents D q of the average generalized inverse 
participation numbers (Z q ). 

In the following section we show that for the octonacci 
sequence and its associated labyrinth tilings the scaling 
exponents D dd of the average generalized inverse partici- 
pation numbers (Z q ) in d dimensions are multiples of the 
one-dimensional scaling exponents D q d . Further, we show 
that for the other quasiperiodic systems this relation is 
approached for large system sizes N a . All following deriva- 
tions can be equivalently applied to the scaling of the av- 
erage participation ratio (p), which is just a function of 
Z 2 , and hence it can be derived that the scaling exponent 
7 is independent of the dimension. 
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( a ) energy E 3 ) energy E 



Fig. 1. Energy dependence of the scaling exponent Di of the inverse participation number Zi(E) averaged over an energy 
interval [E - AE/2, E + AE/2] with AE = 0.05 for the silver-mean chain Cf x g for v = 0.8 (a) and for different values of v (b). 
In (a) also the values of the participation ratio p(^) and the scaling exponent 'y(E) of the averaged participation ratio p(E) as 
well as the density of states (DOS) are shown. 



3.1 Special case - silver-mean model 



labyrinth tiling starting with the general definition of the 
generalized inverse participation numbers: 



We start with the silver-mean model, for which all se- 
quences are palindroms and only even system sizes occur. 
Thus, the lattice resolves into 2 d ~ 1 identical lattices in d 
dimensions [16]. This allows us to decompose the calcula- 
tions of the higher-dimensional inverse participation num- 
bers Z q (<P) for a wave function with arbitrary dimension 
into a product of inverse participation numbers Z q of 
the one-dimensional wave functions W |26| according to: 



z q ($ s ) = E i* 



2<l 



(16a) 
(16b) 



,_, -n N a N a N a 

^^EI^EI^-El^r dec) 



2 d - 



1=1 



n=l 



2tt- l )(i- 1 )z q (& i )Zg{& j ) . . . Z q {& k ) 



(16d) 



(17a) 



sec 



N a /2 N a /2 N a 

E^)---E- 

i=i j—i fe=i 



e ^E^)---E wE^ fe ) 



|3j 2«( rf - 1 ) 1 



(17b) 

N a N a N a 

^E^)---E^(^)E^) 

-1 j — l m—l 

(17c) 

= 2 (d-i)( q -i) {z id )d „ N -d< A _ (y*Ly r F . (i 7d ) 

However, the separation in 2 d_1 identical lattices used 
here cannot be exactly applied to the other sequences. 
In the following we will derive the equations for other 
metallic-mean models. 



The step from Eq. ( 16b I to Eq. ( 16c ) uses the fact that 
the grid is 2 d_1 -partite in d dimensions, i.e. it resolves into 
2^-i identical labyrinth tilings, which only differ by their 
starting point for the construction rule. Consequently, each 
of these tilings contributes with the same amount to the 
inverse participation number. In Eq. ( 16b ) we only sum 



over the vertices of one of th ese tilings. Due to this re- 
lation the summation in Eq. ( 16c ) can be separated into 
2 d_1 parts, which each yield the same result, and hence, 
we obtain the desired product structure. 

Assuming that the average inverse participation num- 
ber obeys the scaling law (Z q d ) ~ N~ T <> with r^ d = 
Dq d (q— 1) for the one-dimensional octonacci chain, we can 
derive scaling expressions for the associated <i-dimensional 



3.2 General case - arbitrary metallic-mean chains 

For some of the metallic-mean models even as well as odd 
chain lengths N occur, where the eigenstates show some 
significant differences for these two models as outlined in 
Sec. [2] Thus, for the proof we distinguish between even 
and odd approximants. 



Even chain lengths N a : We start with the easier case of 
even chain lengths N. In this case we can use an anal- 
ogous derivation as for the silver-mean model. Arbitrary 
metallic-mean chains are no palindroms and, thus, the grid 
in d dimensions does not resolve into identical labyrinth 
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Fig. 2. Scaling exponents 7 of the participation ratio (p) and D2 of {Z2) averaged over the complete energy spectrum in one, 
two and three dimensions for the golden-mean model (a), the silver-mean model (b), and the bronze- mean model (c). 



tilings. However, due to the hidden mirror symmetries of 
the chains as discussed in Sec. l2.2l the wave functions of the 
different labyrinth tilings approach each other for N 3> 1, 
because then the influence of the boundaries and the single 
phason fli p for b > 2 b ecom es negligible. Hence, the step 
from Eq. (16b) to Eq. (16c) is still approximately correct 
for N > 1. 

The deviations of the average participation ratios for 
the two possible tilings in two dimensions, C and £*, are 
shown for different system sizes in Fig. [3j The results are 
quite close for both tilings of the golden- and bronze-mean 
model, respectively. Further, the deviations of the partic- 
ipation ratio for both tilings become smaller for higher 
approximants. To be specific, the total difference of the 
average participation ratio is smaller than 1% for systems 
sizes V4 Au 



and Vo 



= (Nq) 2 /2 = 98 for the golden-mean model 
98 for the bronze-mean model. 



§ 1 



Ni + l 
2 



M-l N a N a N a 

£]>> 9 (^) + E E 

i—l 3 — 1 i=M+l 3=1 

M-l N a 

E^(* Mi )+ E Z ^ M3 ) + 

3=1 j=M+i 
Z q ($ MM ) 



+ 



{p L_ 

Ni + 



- V V v l#? \ 2q + 2 V \$f* M 

1 Z^ Z^ Z^ I Im | N 2 + 1 Z-^ 



(18c) 

2q 



i=l 3=1 l.mGC 



Na>l 1 



1 E E 



2 E 1*1 



JV2 + 1 



E KT + E k 



E |*J 

|2« 



(18d) 

,2q 



+ 



MM 
Im 



(18e) 



Odd chain lengths N a : In the next step we have to derive 
such a relation also for odd system sizes. In this case the 
proof is more complicated because we have to treat the 
state E M = separately due to the special structure of the 
corresponding wave function (cp. Eq. Q). In the following 
we show the derivation for the two-dimensional case, for 
which the grid decomposes into the two tilings C and £*, 
where the state E M — belongs to the first one. This 
reasoning can be extended to arbitrary dimensions in a 
straightforward way. 

Starting with the definition of the average generalized 
inverse participation number we obtain the scaling law 
from the following derivations: 



N a N a N a N a 



"M-l N a M 

E E z «(* iJ ) + E z i( $Mj ) 

. i=l 3=1 3 = 1 



Vi 

2 

N2 + 1 



(18a) 
(18b) 



a i=i j = i 1=1 m=l 



N a N a 



2 8-l 



N2 + 



;=i m=l 

N a N a N a N a 

T EEEE|^r+ 



i=l j = l 1 = 1 m=l 
N a N a 



o 09— 1 

TV J + 1 ^ ^\ 1 " 



2,, 



(=1 m =l 

N a N a 



(18f) 



(18g) 



' 1 1 3 = 1 m=l 



8=1 !=1 
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(18h) 
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Fig. 3. Difference of the participation ratios p(<P a ) for the two 
possible labyrinth tilings, C and C , for the golden- and bronze- 
mean model in two dimensions for a coupling parameter v = 
0.3. Results are shown for C% u with N = 14 (a), £$ u with 
N = 35 (b), £f r with N = 14 (c), and £? r with N = 44 (d). 



N a 



(v a 2d y rlq 



(18j) 



The step from Eq. (18b I to Eq. (18c) is possible due 
to the same argu ment s as in the case with even chain 



lengths N a . In Eq. ( 18f I we have to take into account that 
the state i = j = M only differs from zero either on C 
or £*. In our case this means that only the product of 
the state of Eq. Q is considered for the state E MM 
because W M ~ vanishes on the labyrinth tiling C. All other 
steps of the proof follow the same way as for even system 
size N a . Analogously the same concept can be applied to 
higher dimensions. 

Hence, we have shown that the scaling of the average 
inverse participation numbers (Z q ) with the linear system 
size N 3> 1 is described by the one-dimensional scaling 
exponent according to 



DT = dD , 



Id 



(19) 



Since the system size grows also with Vdd ~ N d , the eigen- 
states in these metallic-mean systems maintain their frac- 
tal nature, i.e. wave functions possess the same scaling 
exponent in arbitrary dimensions with respect to the ac- 
tual system size V^d- 



4 Conclusion 



In this paper we have investigated the spatial distribution 
of the wave functions for the generalized labyrinth tilings 
in different dimensions d. First we have introduced the 
construction rules for the eigenstates in labyrinth tilings 
in two and three dimensions. The numerical results for 
the scaling behavior of the generalized inverse participa- 
tion numbers Z q revealed the multifractal properties of 
the wave functions. In particular, we found numerical ev- 
idence that the spatial distribution of the wave functions, 
which is described by the scaling exponents 7 and D q , 
is independent of the dimension d. This result has been 
substantiated by showing analytically that for the silver- 
mean model the scaling exponents D q in d dimensions 
are always multiples of the one-dimensional scaling expo- 
nents D q d according to Eq. (19 1. For the other metallic- 



mean chains we have shown that the scaling exponents 
approach this relation asymptotically. This result shows 
that the wave functions in the metallic-mean systems are 
characterized by the same spatial expansion and the same 
multifractal properties in each dimension. Hence, the gen- 
eralized labyrinth tilings with their special structures can 
be useful toy models to study the characteristics of higher- 
dimensional quasiperiodic systems efficiently by numerical 
methods and also analytically. 
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